Anomalous magnetotransport in (Yi_ a; Gd 2: )Co2 alloys: interplay of disorder and 

itinerant metamagnetism. 

A. T. Burkov, A. Yu. Zyuzin 
A. F. Ioffe Physico- Technical Institute, Russian Academy of Sciences, Saint- Petersburg 194021, Russia. 



o 
o 

(N 
S3 



c3 



T3 
S3 

O 

i o i 

(N 
> 
On 

m 
in 

o 
m 
O 

03 



i 1 

■ 

c : 
o 

o ■ 



X 



T. Nakama, K. Yagasaki 
Physics Department, University of the Ryukyus, Okinawa 903-0213, Japan 
(Dated: February 2, 2008) 

New mechanism of magneto-resistivity in itinerant metamagnets with a structural disorder is in- 
troduced basing on analysis of experimental results on magnetoresistivity, susceptibility, and mag- 
netization of structurally disordered alloys {Y\- x GA x )Co2- In this series, YC02 is an enhanced 
Pauli paramagnet, whereas GdCo2 is a ferrimagnet (T c =400 K) with Gd sublattice coupled antifer- 
romagnetically to the itinerant Co-3d electrons. The alloys are paramagnetic for x < 0.12. Large 
positive magnetoresistivity has been observed in the alloys with magnetic ground state at temper- 
atures T<T C . We show that this unusual feature is linked to a combination of structural disorder 
and metamagnetic instability of itinerant Co-3d electrons. This new mechanism of the magnetore- 
sistivity is common for a broad class of materials featuring a static magnetic disorder and itinerant 
metamagnetism. 

PACS numbers: 71.27.+a, 72.15.Gd, 75.10.Lp, 75.30.Kz 



I. INTRODUCTION 

Interplay of structural disorder and magnetic inter- 
actions opens a rich field of new physical phenom- 
ena. Among them are the actively discussed possi- 
bility of disorder-induced Non-Fermi Liquid (NFL) be- 
havior near a magnetic Quantum Critical Point (QCP) 
as well as a broader scope of effects of disorder 
on magnetotransportiSii Structurally disordered alloys 
Yi_2;Gda;Co2 are quasi-binary solid solutions of Laves 
phase compounds YC02 and GdCo2 . The compounds be- 
long to a large family of isostructural composites RC02. 
YC02 is an enhanced Pauli paramagnet whose itinerant 
Co-3d electron system is close to magnetic instability. In 
external magnetic field of about 70 T this system under- 
goes a metamagnetic transition into ferromagnetic (FM) 
ground stated GdCo2 is, on the other hand, a ferrimag- 
net with a Curie temperature of 400 K in which the 
spontaneous magnetization of 4f moments is anti-parallel 
to the induced magnetization of the Co-3d band. Com- 
pounds of this family and their alloys provide a conve- 
nient ground for experimental studies of magnetotrans- 
port phenomena. The electronic structure in the impor- 
tant for the transport vicinity of the Fermi energy is com- 
posed mainly of Co-3d states and is, to the first approx- 
imation, the same for all compounds of RC02 family. It 
has been found that the main contribution to the resis- 
tivity of RC02 comes from the scattering of conduction 
electrons on magnetic fluctuations due to strong s-d ex- 
change coupling,^ therefore the transport properties are 
expected to be especially sensitive to the magnetic state 
of the sample. GdCo 2 occupies a special place in RC02 
family since Gd 4f magnetic moment has no orbital con- 
tribution and, therefore crystal-field effects are not im- 
portant for this compound. 

The experimental results on the transport properties 



of Yi_ 2; Gd : ];Co2 alloys has been published partly in our 
previous articled Here we analyze these and new experi- 
mental results in order to reveal the physical mechanism 
of anomalous megnetotransport properties observed in 
the alloys. In this paper we will discuss the magneto- 
transport properties of the FM alloys. The properties of 
paramagnetic alloys will be published elsewhere. 



II. EXPERIMENTAL 

Samples of Yi_. r Gd a ;Co2 were prepared from pure 
components by melting in an arc furnace under a protec- 
tive Ar atmosphere and were subsequently annealed in 
vacuum at 1100 K for about one week. An X-ray analy- 
sis showed no traces of impurity phases. A four-probe dc 
method was used for electrical resistivity measurements. 
Magnetoresistivity (MR) was measured with longitudi- 
nal orientation of electrical current with respect to the 
magnetic field. The size of the samples was typically 
about lx 1x10 mm 3 . Magnetization was measured by 
a SQUID magnetometer for samples from the same in- 
got as that used for the resistivity and AC susceptibility 
measurements. 



III. EXPERIMENTAL RESULTS 

The magnetic phase diagram of the Yi_ a; Gd a ;Co2 
system inferred from the transport and magnetic 
measurements 6 is shown in Fig. ^ Curie temperature 
T c decreases with increasing content of Y and eventu- 
ally drops to zero. A precise determination of the crit- 
ical concentration x c which separates the magnetically 
ordered ground state and the paramagnetic region is dif- 
ficult, since on the onset of the long range order its signa- 
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FIG. 1: The upper panel shows the ordering temperature 
T c ■ (right y-axis), and the MR (left axis) of the 
Yi-zGdzC^ system*. The MR was measured at T = 2 K 
in magnetic field of 15 T. The dotted vertical lines indicate 
phase boundaries at zero temperature. The lower panel dis- 
plays normalized resistivity 
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tuxes in the magnetic and transport properties are very 
weak. The first firm evidence of the long range order are 
found for alloy with x=0.14 in ac susceptibility at T=27 
K, Fig. El 

Quantum critical scaling theory predicts that when 
Curie temperature of a FM system continuously depends 
on an external parameter x, this dependence is expressed 

2 

as*: 



T c oc \x - x c \ d +'- 2 

with critical index z = 3 for a FM system of spatial 
dimension d = 3. The experimental T c vs. x depen- 
dency does follow this relation, but with additional kink 
at x = Xt . A possible origin of this kink will be discussed 
later. Linear extrapolation of the phase separation line 
on the phase diagram Fig. 2] to T c = gives as the crit- 
ical concentration x c — 0.12. We do not claim however 
that QCP exists in this alloy system. Direct experimental 
verification that T c — > as x approaches x c from mag- 
netically ordered state is difficult for a disordered alloy 
system. 
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FIG. 2: The ac susceptibility of the Yi„ :E Gd :E Co2 alloys. 
Note, the experimental data for YC02 and for the alloy with 
x — 0.14 are multiplied by factor 20. 



The very surprising result is the positive MR in the 
FM phase at low temperatures, see Figs. ^ and [3J The 
well known theoretical result for MR of a localized mo- 
ment ferromagnet was derived long ago by Kasuya and 
De Gennesi As it follows from their theory, MR of a 
metallic ferromagnet should be negative, having a maxi- 
mum absolute value at Curie temperature, and approach- 
ing zero as T — > 0, and in the limit of high temperatures. 
Qualitatively this behavior has been supported by ex- 
periment, as well as by later more detailed theoretical 
calculations. The present experimental results are in a 
qualitative agreement with this theoretical behavior only 
for alloys with large Gd content (x ^ 0.4) (Fig.©. MR 
of the FM alloys with composition 0.3 > x > 0.14 fun- 
damentally differs from this theoretical behavior. Let us 
note that this composition range falls into the region of 
the phase diagram between the paramagnetic phase and 
the additional phase boundary indicated by the kink in 
the T c vs x dependency, see Fig. ^ MR of these alloys 
is positive below Curie temperature and is very large. 

The known mechanisms of a positive MR can not ex- 
plain the experimental data. A rough estimate of Lorenz 
force-driven MR one can get from a comparison with the 
MR of pure YCo2^ In the most pure samples of YC02 
(with residual resistivity of about 2 /if2cm) the Lorenz 
force-driven positive contribution to the total MR does 
not exceed 5%. On the other hand the resistivity of the 
FM alloys at low temperatures falls into the region from 
30 to 100 /xS7 cm, i.e. at least one order of magnitude 
larger than the resistivity of pure YC02. Therefore, ac- 
cording to Koehler's rule, this mechanism can give MR 
of only about 0.5%, this has to be compared with the 
experimental MR of almost 40%. 

Weak localization effect is known to give positive MR. 
However our estimates show that this mechanism can 
give a contribution which is at least two order of magni- 
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spatially fluctuating effective field induces an inhomoge- 
neous magnetization of Co-3d electron system: 
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FIG. 3: The upper panel shows the temperature dependence 
of the MR of the Yi_ 2 ,Gd I Co2 alloys, measured in field of 
15 T. Large positive MR of FM alloys (x < 0.3) is observed 
at low temperatures. The field dependencies of MR, measured 
at T=2 K, are presented on the lower panel. 



tude smaller than the observed MR. 



IV. DISCUSSION 

The key for understanding the mechanism of the pos- 
itive MR is a combination of strong dependence of the 
magnetic susceptibility x of metamagnetic Co-3d subsys- 
tem on the effective magnetic field and of the structural 
disorder in the R - sublattice of the alloys. In case of 
GdCo2 as well as in the case of other magnetic RC02 
compounds with heavy R-elements, the 4f-3d exchange 
interaction is described by introducing an effective field 
which acts on 3d electrons as: 



B, 



eff 



rifd-Mf — B 



where B is external field, My is the uniform magneti- 
zation of R - sublattice, and rifd - is the f-d coupling 
constant (in case of GdCo2 rtfd ~ 50 T/f.u./zg^). In 
Yi_ 2; Gd 2 ;Co2 alloys the Gd moments are randomly dis- 
tributed over the R-sites of the crystal lattice. Therefore, 
the effective field acting on 3d electrons depends on the 
local distribution of Gd moments and is therefore a ran- 
dom function of coordinate. This random field can be 
characterized by a distribution function P{_B e g(r)}. The 



m(r) = x (Bcs) B cS (r) . 



Therefore even at zero temperature in the ferromagnetic 
ground state, there are two kind of static magnetic fluc- 
tuations in the system: i. Mf(r); ii. m(r). These fluctu- 
ations give an additional contribution to the resistivity. 
Since at T=0 K in ferromagnetic phase the 4f magnetic 
moments are saturated the corresponding contribution 
to the resistivity does not depend on external magnetic 
field. On the other hand, the 3d magnetic moment, as 
we will see later, is not saturated even in the ferromag- 
netic ground state. Therefore the 3d magnetization does 
depend on the external field. However, as long as 3d sus- 
ceptibility is field independent and uniform, the external 
magnetic field will change only the mean (non-fluctuating 
part) value of the magnetization, whereas the magnitude 
of the fluctuations of m and corresponding contribution 
to the resistivity remain unchanged. However actually, 
the 3d system is close to the metamagnetic instability. 
Therefore the 3d susceptibility x is field dependent and 
gives rise to static magnetic fluctuations which are de- 
pendent on magnetic field resulting in non-zero magne- 
toresistivity. 

For a qualitative analysis of this new mechanism of 
MR we make the following assumptions. 

1. We consider the system only in its ground state. 

2. Correlations between potential and spin-dependent 
scattering are neglected. 

3. Only the contributions which may be strongly depen- 
dent on external magnetic field are retained, e.g. we do 
not include the effects related to a change of d-density 
of states in magnetic field, and contributions due to 
potential scattering and scattering on 4f moments. In 
fact only the contribution related to the scattering on 
fluctuations of Co-3d magnetization will be considered. 

Figure 0] shows schematically the dependency of the 
Co-3d magnetization of RC02 compounds on effective 
magnetic field 9 and the distribution function P{B c s(r)}. 
The metamagnetic transition is indicated by the rapid in- 
crease of the magnetization around 70 T. For the further 
discussion it is important to have an estimate of the mag- 
nitude of the fluctuations of the effective field, i.e. the 
width of the distribution function P{B c g(r)}. We can 
get a hint considering experimental data on field (0 to 
7 T) dependency of the magnetization, FigJSJ 

Two important points follow from these data: i. there 
is a non-saturating para-process in these field dependen- 
cies above about 2 T; ii. the estimated susceptibility of 
this para-process (« 0.016 //b/T) is larger than the sus- 
ceptibility of 3d system below and above metamagnetic 
transition (« 0.002 /xs/T), however it is smaller than the 
susceptibility in the transition region (« 0.04 /xg/T)£ 
This indicates, that the scale of the fluctuations is larger 
than the width of the metamagnetic transition. There- 
fore we can not treat the fluctuating part of the effec- 
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FIG. 4: Dependence of the magnetization of Co-subsystem 
in RC02 compounds vs effective magnetic field (left axis)- 
(shown schematically). The dotted line is a schematic repre- 
sentation of the distribution function P(B e g) (right axis), the 
shaded area indicates the position of P(B) in external field of 
15 T. 



becomes larger than Bq, i.e. at this composition almost 
whole volume is occupied by the high magnetization com- 
ponent, i.e. y as 1. With a further increase of Gd con- 
tent the mean magnetization of Co should increase with a 
smaller rate, determined by the slope of the m(B) depen- 
dency above Bq, however y = 1 in this region. According 
to this scenario p m will increase with x at first, reach a 
maximum value at x which corresponds to y = 0.5 and 
will decrease with further increase of x approaching to 
zero at x ss x t . We believe that x t corresponds to the 
kink on the phase diagram, Fig. ^ The expected vari- 
ation of pm with the alloy composition is schematically 
shown in Fig. The total experimental resistivity in- 
cludes additionally contributions coming from potential 
scattering and from scattering on 4f moments (po)i which 
both are proportional to x(l — x) and also are depicted 
in Fig. For comparison the experimental resistivity is 
shown in this picture too. We present here the resistiv- 
ity measured at T=2 K in field of 15 T to exclude the 
contributions related to spin-flip scattering. 




FIG. 5: The magnetization of the Yi_ a; Gd :E Co2 alloys vs mag- 
netic field. 

tive field as a perturbation and have to resort to a phe- 
nomenological analysis. 

The distribution function P(B e s) describing the fluc- 
tuating effective field depends on the alloy composition. 
For diluted alloys (x ss 0) the most probable value B av of 
B e g is close to zero. As x increases B av shifts to higher- 
values and in a certain range of the concentrations, the 
function will have essentially non-zero weight for both 
£> e ff < Bq and £? e ff > Bq, see Fig. In this case there 
shall be regions with low and with high 3d magnetization 
in the sample. The resistivity resulting from this static 
disorder in 3d magnetization can be expressed as: 

p m = p S d-y{i-y) ■ 

Where y = J^ g P (B e g) dB is the volume fraction of the 
high magnetization component. Parameter y depends on 
the alloy composition x, and on external magnetic field 
B. In zero field, point y — corresponds to x = 0. 
As the content of Gd increases B av shifts to larger effec- 
tive fields, and finally, at some alloy composition x tl it 
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FIG. 6: Schematic of contributions to the low-temperature 
resistivity of Yi_ a; Gd :!; Co2 alloys. • -experimental normalized 
resistivity measured at T=2 K in field of 15 T. Solid line 
represents schematically p m , broken line - po. 

The experimental low-temperature resistivity shows 
the expected variation with x (one needs to keep in 
mind that the relation between x and y is in general 
non-linear, especially around magnetic phase boundary 
x c ). The resistivity attains the maximum value of about 
100 /if2 cm in the region, which corresponds to maximum 
static magnetic disorder at y « 0.5 (room temperature 
resistivity of the alloys weakly depends on x and is about 
150 pVlcm). About the same value was obtained for the 
high temperature limit (maximum magnetic disorder) of 
magnetic part of the resistivity, arising from scattering on 
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3d temperature- induced magnetic fluctuations in YCo2^ 
This suggests that the main part of the experimental 
low temperature resistivity at x < x t originates from the 
scattering on the magnetic fluctuations of 3dmagnetiza- 
tion, i.e. is identical to p m . 

Let us show that p m is of order of the spin-fluctuation 
contribution to resistivity at large temperatures, i.e. is of 
order of 100 /ifkm. In case of YC02 the spin-fluctuation 
contribution is the most important and is nearly indepen- 
dent of temperature above about 200 K. Hamiltonian of 
s-d exchange interaction is given by 



H sc1 = gJ drs(r)S d (r) 



here s(r) and Sd(r) are spin density of s- and d- elec- 
trons, correspondingly^ 

Spin fluctuation contribution to resistivity has formifi: 



ttG 2 N s — / dqq 3 / n 

4ne 2 S TJ HH J sinh 2 (w/2T) 

—00 



_ / 2k F 
X \ 9 ~kf' W 

where N s is density of states of s-electrons, |f- is ratio 
of Fermi momentum of s- and d- electrons. Dynamic 
susceptibility x (<?) w ) is given by the equation: 



x" 1 ^) = x ~ 1 («)(i -Wr,) 

here T q is damping of the spin-fluctuations, whereas 
static nonlocal susceptibility x 1 (<?) is given by 



X- 1 (q) = X' 1 + Aq 2 /N d 

with x — Nd/((T) , where C(^) is inverse Stoner fac- 
tor, renormalized by spin fluctuations, Nj is density of 
states of d-electrons, and A < 1 is a dimensionless con- 
stant. 

At large temperatures T > T q the expression Q re- 
duces to: 



3nm 
Ane 2 



G 2 N S T X 



(2) 

The last equality is valid when 2Ak F lk* F < ((T) < 1 
which must be the case for YC02. Note that neglecting 
momentum dependence of susceptibility means that s- 
electrons see the d-spin fluctuations as point scatterers. 

Scattering of conducting s - electrons by the static 
random distribution of spin density of d-electrons we can 
estimate in the following way. 

Due to random distribution of magnetic moments of 
Gd, s-electrons experience scattering by 

j (Sd <?)) = - J Xm (r-T')SB eff (r') 



random potential. Here Xm( r — r ') is nonlocal suscepti- 
bility near metamagnetic transition. We estimate corre- 
lation function of fluctuating effective field as 

{SB eff (r)SB eff (r')) = (2S Gd n id ) 2 S (r - r') x (1 - x) a 3 . 

Here a 3 is volume of the formula unit. In Born approxi- 
mation the corresponding contribution to the resistivity 
is given by the expression: 



1 

Pm = ^G 2 N S (2S Gd n id ) 2 x (1 - x) a 3 f dqq 3 X 2 m (q^f 
ne J \ ftp 



(3) 

Both expressions {Q and © give the contributions to 
the resistivity due to scattering on spin fluctuations, how- 
ever in the first case they are of thermal origin, whereas 
in the second case the fluctuations are due to randomness 
of the effective field. 

Assuming that nonlocality of the susceptibility is not 
important we obtain from J5J and (0: 



p m (2S Gd n id ) x (1 - x) a 3 x 2 m 



p 3ttTx 

Using experimental results for Xm and T^iaLii we find 
^ in the range from 0.5 to 3, i.e. the resistivity caused 
by the static magnetic fluctuations is of the same order 
as the temperature-induced spin fluctuation resistivity. 
The uncertainty is mainly due to determination of Xm 
near the metamagnetic transition. Taking Xm as the sus- 
ceptibility of YC02 at B=70 T (at the field of metamag- 
netic transition)^ gives the upper bound for Bm. , Whereas 
Xm for the disordered alloy of x=0.18, estimated from our 
results on M(B), Fig. [S] gives the lower bound. 

In external magnetic field, the effective field _B e ff 
decreases^, therefore y also decreases. Depending on 
the value of yo - the volume fraction in zero field, p m 
will cither increase or decrease, resulting in positive or 
negative magnetoresistivity: for 0.5 < yo < 1 it will 
be positive, whereas for < yo < 0.5 we will have a 
negative MR. In agreement with the model, the exper- 
imental MR is positive at 0.15 < x < x t and quickly 
decreases at x > Xt = 0.3 where y ~ 1. Nearly lin- 
ear field dependencies of MR, observed for x < 0.3, see 
Fig. |3 implies that the width of P(B e g) in this compo- 
sition range is larger than our experimental field limit 
of 15 T. The region y < 0.5 almost coincides with the 
paramagnetic region of the phase diagram. The model 
predicts negative MR for this region, and this predic- 
tion agrees with the experimental result. The model also 
gives a satisfactory description of the residual resistivity 
behaviour in this region, see Fig- EI This agreement sug- 
gests that at x < x c the system is actually in spin-glass 
state. In this region additional essential contributions to 
MR are present. First, in the paramagnetic region there 
is a negative MR due to suppression of magnetic disor- 
der in 4f magnetic moment system by external magnetic 
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field. This negative contribution may be the reason why 
the cross-over point from the positive to negative mag- 
netoresistivity does not coincide with the maximum of 
resistivity, see Fig. ^ Secondly, there can be additional 
contributions, both positive and negative, near to zero- 
temperature magnetic phase boundary due to closeness 
to QCP. 

An independent test of the model is based on the ob- 
servation that the critical magnetic field Bq of the meta- 
magnetic transition in 3d subsystem increases under ex- 
ternal hydrostatic pressure (P)>i2*i2i Therefore, basing on 
the model, we expect that the resistivity of the alloys 
with the composition left of the resistivity maximum (see 
Fig.nj will decrease with pressure, whereas for the alloys 
right of the maximum it will increase with the increas- 
ing pressure. The experimental results for three alloy 
compositions are shown in Fig. [7| The sign of the pres- 




Pressure, kBar 

FIG. 7: The resistivity of the Yi-^Gd^Cc^ alloys vs external 
hydrostatic pressure at T = 2 K. 



sure effect is in agreement with the model prediction: the 
resistivity decreases with pressure for x=0.l, whereas it 
increases with P for £=0.18 and x=0.3. Moreover, there 
is a good scaling of pressure and magnetic field depen- 
dencies of the resistivity, Figs. |21 and (P— > aB) with 
scaling parameter a which is close to literature data on 
the pressure dependence of B : ^fr ~ 1.5 T/kBarii2i±£ 



V. CONCLUSION 

It has been found that at low temperatures there is 
a large contribution to the resistivity related to scat- 
tering on magnetic fluctuations in metamagnetic itin- 
erant 3d system, induced by fluctuating effective field 
of 4f moments. Large positive MR, found in the FM 
Yi_ a; Gd 2 ;Co2 alloys and strong pressure dependence of 
the resistivity are explained as arising from a combina- 
tion of static magnetic disorder and strong magnetic field 
dependence of magnetic susceptibility. 

We want to emphasize that this mechanism of resistiv- 
ity (and of MR) is not material specific, rather it should 
be common for a broad class of disordered itinerant meta- 
magnets with strong coupling of conduction electrons to 
the magnetic fluctuations. In Yi_ a; Gd a ;Co2 the relevant 
disorder originates from random distribution of d-f ex- 
change fields, however a similar effect should arise when 
there is a random distribution of local susceptibilities (a 
corresponding treatment for Kondo systems was recently 
developed wM>). Positive MR observed in FM alloys 
Y(Coi_j; A\ x )^^ and in Er^Yi-^Coai^ may be explained 
by this mechanism. 
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The effective field decreases for antiferromagnetic 4f-3d ex- 
change, for ferromagnetic exchange (like that in RC02 com- 
pounds with light R-elements) the effective field increases 
with the external field. 



